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ABSTRACT 

We examine the effects of low-temperature, or An = 0, dielectronic 
recombination (DR) on the ionization balance of the Fe M-Shell (Fe IX through 
Fe XVI). Since An = rates are not available for these ions, we have derived 
estimates based on the existing rates for the first four ionization states of the 
CNO sequence and newly calculated rates for L-shell ions of 3rd row elements 
and Fe. For a range of ionization parameter and column density applicable 
to the intrinsic absorbers detected in ASCA, Chandra, and XMM-Newton 
observations of Seyfert galaxies, we generated two grids of photoionization 
models, with and without DR. The results show that the ionization parameter 
at which the population of an Fe M-shell ion peaks can increase in some cases 
by factor > 2 when these rates are included. More importantly, there are 
dramatic changes in the range in ionization parameter over which individual 
M-shell ions contain significant fractions of the total Fe (e.g. > 10%) in the 
plasma. These results may explain the mismatch between the range of Fe 
ionization states detected in the X-ray spectra of Seyferts, identified by the 
energies of the M-Shell Unresolved Transition Array, and those predicted by 
photoionization models of the X-ray absorbers that reproduce lines of second 
and third row elements. The results suggest that care should be taken in using 
3rd and 4th row ions to constrain the physical conditions in photoionized X-ray 
plasmas until accurate DR rates are available. This underscores the importance 
of atomic physics in interpreting astronomical spectroscopy. 
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1. Introduction 

New X-ray observatories such as the Chandra X-ray Observatory and XMM- Newton 
have made high resolution spectroscopy of active galaxies possible for the first time, allowing 
the detection of phases of gas that were previously unknown. The physical state of this 
gas is interpreted by reference to large-scale numerical simulations, hence, it is critical that 
physical processes that occur in an ionized plasma be fully understood (the main topic of 
the conference "Spectroscopic Challenges of Photoionized Plasmas" , held at the University 
of Kentucky, Lexington, Kentucky, USA on the 15 - 18 November, 2000; [Ferland & Savin 
2001]). 

The ionization structure and temperatures in cosmic plasmas depend strongly 
on ionization and recombination processes. Recombination can occur via radiative 
recombination, charge transfer, or dielectronic recombination (DR). The DR process 
involves the capture of a free electron by an ion, resulting in doubly-excited recombined 
species with the electron in an auto-ionizing state. This is followed either by auto-ionization 
(i.e., the reverse process) or the emission of a photon, the latter leaving the recombined ion 
in a stable state. 

There are two varieties of DR. The first involves the capture of the electron to 
auto-ionizing states with high principal quantum number, with an accompanying core 
excitation with a change in principal quantum number An > (e.g. Burgess 1964). This 
form of DR dominates over radiative recombination at temperatures comparable with 
the ionization potential, i.e coUisionally ionized plasmas. Nevertheless, terms of the same 
parity as the ground term exist between the ground state and the first excited state of 
the opposite parity, and thus a series of autoionizing states converge on these terms (see 
Nussbaumer & Storey 1983). This results in a second type of DR in which there is no 
change in the principal quantum number of the excited core, or An = 0, which can proceed 
at low temperatures for which the thermal energy of the free electron is less than the 
excitation energy of first resonance of the recombining ion. Low-temperature DR is usually 
the dominant process in photoionized plasmas, with temperatures typically 5 - 10% of 
the ionization potential, such as those detected in X-ray spectra of Seyfert galaxies (e.g 
Reynolds 1997; George et al. 1998; Sako et al. 2000; Kaspi et al. 2001, 2002; Kinkhabwala 
et al. 2002; Turner et al. 2003). 

While low-temperature DR rates have been calculated for a number of abundant 2nd 
and a few 3rd row elements (Nussbaumer & Storey 1983; Nussbaumer & Storey 1986; 
Nussbaumer & Storey 1987), the full set of ionization states have only been completed for 
2nd row elements. One of the difficulties in these calculations is that the rates depend on 
the accurate determination of the low-level autoionizing states (see Gu 2003). Ali et al. 
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(1991) estimated low-temperature DR rate coefficients for other ions by taking the mean 
values from the first four ions of C, and O calcTilatcd by Nussbaumer & Storey (1983), 
and demonstrated that the model-predicted emission-line strengths were strongly affected 
when these rates were introduced. Savin (2000) has shown that uncertainties in DR rates 
translate into similar uncertainties in predicted ionic column densities, which may hamper 
the determination of elemental abundances ratios in the intergalactic medium. There has 
been recent work in determining both high- and low- temperature DR rates for L-shell (the 
Li- through Ne-like iso-electronic sequences), both experimentally (Savin et al. 1999, 2002) 
and computationally (Gu 2003; Badnell et al. 2003), the results of which emphasize the 
importance of this process in photoionized plasma. Note that, due to the difficulties in 
calculating the energies of the autoionizing states, experimental results are critical. 

On the other hand, there has not been any recent work published on An = rates 
for M-shell ions (the Na- to Ar-hke sequence). This is unfortuitous, since M-shell Fe ions, 
Fe IX to Fe XVI, contribute strongly to the unresolved transition array (UTA) of inner shell 
n = 2 — 3 lines seen in absorption in high resolution X-ray spectra of several Scyfcrt galaxies 
observed with XMM-Newton and Chandra such as IRAS 13349+2438 (Sako et al. 2001), 
NGC 3783 (Blustin et al. 2002; Kaspi et al. 2002; Netzer et al. 2003), Mrk 509 (Pounds et 
al. 2001; Yaqoob et al. 2003), and NGC 5548 (Steenbrugge et al. 2003; Steenbrugge et al., 
in preparation). For N- and M-shell Fe, the UTA features he between ~ 15.5 A and 17 A , 
depending on which ions are present (Behar, Sako, & Kahn 2001). In principal, the UTA 
can be used to constrain the ionization state of the absorbers, which is especially valuable 
if the resonance lines from H-, He-, and Li-hke C, N, and O are saturated. However, in 
the case of NGC 3783 (Netzer et al. 2003), the ionization state suggested by the UTA 
wavelengths seems lower than that predicted by photoionization models generated to fit 
CNO absorption features and the ionization state/column density of 3rd row elements such 
as Si and S. A similar discrepancy appears to be the case for NGC 5548 (Steenbrugge et al., 
in preparation). One of the possible causes Netzer et al. suggested to explain this mismatch 
is the lack of An = rates for the M-shell Fe ions, which leads to the models predicting the 
Fe in too high an ionization state (compared to lower Z elements). 

In order to explore the effect of An = DR on the M-shell Fe ionization balance, we 
adopted the approach described in Ali et al. (1991), adding estimated rate coefficients to 
the development version of the photoionization code Cloudy (Ferland et al. 1998). In order 
to gauge the magnitude of the effect, model grids were generated with and without these 
rates. We also examined the effects of uncertainties in the estimates rates. The model 
details and results are described below. 
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2. Modeling Methodology 

Since the mean An = rate coefficients for CNO show a obvious dependence on the 
residual charge of the ion before recombination (Ah et ah 1991), we would expect a strong 
charge dependence among the 3rd and 4th row M-shell ions. In Figure 1, we show the 
rates for the L-Shell Fe and Ar ions from Gu (2003) for an electron temperature of 4 eV 
(~ 5 X lO^K), similar to what we would expect in photoionized plasmas in which the iron 
abundance was peaked at the M-shell. Since they span the same range in residual charge, 
one might expect that the M-shell Fe rates would be similar to the L-shell Ar rates. Hence, 
we estimated the rate coefficients for iron by generating 2'^'^-order polynomial ffis to the 
L-shells rates of Fe and Ar and the mean rates for the CNO sequence (Ali et al. 1991), with 
the functional form: 

where Z is the residual nuclear charge of the ion before recombination. The individual ffis 
are shown in Figure 1, and the values of the ffi coefficients are listed in Table 1. Since Gu's 
calculations show a different dependence on Z for different elements within the L-shell, 
we opted not to use this fit for the other 3rd and 4th row elements, (in some cases, the 
behavior within the shell cannot be fit with a second order polynomial). However, for the 
sake of completeness, we did estimate An = rates for these ions using a linear fit to set of 
data including the Fe and Ar L-shell rates and the mean rates for the CNO sequence (Ali 
et al. 1991), as shown by the dotted line in Figure 1, with the functional form: 

i? = 3.0 X 10"^^ 10°-^^ cm^ s-\ (2) 

Since An = DR are not fast for recombination from closed shells over the range of 

temperatures we consider, we did not include estimated DR for those ions. Note that we 
did not include temperature dependence for our estimated rate coefficients. 

To quantify the effects of these missing DR rates we performed a series of photoionization 
calculations using the best existing data, but with and without these estimates of the 
An = DR. We did not introduce new high-temperature DR rates, using, instead, the ffi 
to the rates for Fe given in Arnaud & Raymond (1992). For the models, we assumed a slab 
geometry, solar abundances (e.g. Grevesse & Anders 1989), a number density ylh = 10^ 
cm~^ and total hydrogen column density of Nh = 1 x 20^^ cm~^, typical of X-ray absorbers; 
note that the predicted ionization structure is insensitive to uh, unless the densities exceed 
10^° cm~^. For the spectral energy distribution of the incident continuum radiation, we 
assumed a power law of the form F oc where a — 1 ior hv < 13.6 eV, 1.4 for 13.6 eV 
< hu < 600 eV, and 0.8 for hv > 600 eV. We generated the two sets of models over a range 
in ionization parameter (number of ionizing photons per nucleon at the illuminated face 
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of the slab) log U = —1.0 to 1.0, in increments of 0.05 dex. Although the mean electron 
temperatures range from 1.8 x lO^K to 4.5 x 10^ K, the M-shell Fe ions peak between log U 
= —0.5 and 0.5, or 3 x lO^K ^ T ^ 9.0 x 10^ K, justifying the energy we used in deriving 
the rate coefficients. Over this temperature range, the DR rates may vary by factors of ^ 
a few, which is within the uncertainty quoted above. Note that, while T increased slightly 
when the estimated DR rates are included, the ionic columns for 1st and 2nd row elements 
remained essentially unchanged. 

The predicted ionic column densities for Fe Vlll through Fc XVI are shown in Figures 
2 through 4. Since intrinsic X-ray absorption is typically characterized by strong oxygen 
bound-free edges and absorption lines, we also include the O VI, O VII, and O VIII columns 
for comparison. As shown in Figure 2, the ionization parameters at which the column 
densities of Fe VIII, Fe IX, and Fe X peak have all shifted to higher values of U. In the 
case of Fe X, U has shifted by a factor of ~ 1.5. More importantly, the range in ionization 
parameter over which a significant fraction of the total Fe exists as Fe X (i.e., 10%, 
corresponding to a column density ^ 10^^'^ cm~^) has changed from —0.8 < logU < —0.3 to 
—0.8 < logU < 0.1. This suggests that, due to An = DR, relatively low ionization states 
of Fe can contribute significantly to opacity of more highly ionized plasmas than would 
be predicted using the currently available atomic data. Although the peak shifts are not 
significant for Fe VIII and Fe IX, the curves show substantially more of these ions at higher 
values of U. Note that, since Fe IX has a closed shell (Ar-hke), there is no low temperature 
An = process for recombination to Fe VIII, which is produced primarily by radiative 
recombination. The increase in Fe VIII at higher U is due to the larger fraction of Fe IX, 
which is produced by DR from Fe X. The structure in the high ionization parameter tails 
of Fe IX and Fe X is the result of high temperature DR, although this effect is swamped for 
Fe X in the models in which An = DR inchidcd. Fe XI through Fe X III are shown in 
Figure 3, and for each ion, the peaks shift to higher ionization by factors of > 2, reflecting 
the large An = rates. Fe XIV through Fe XVI are shown in Figure 4, and the inclusion 
of An = DR in this case decreases the fractional abundances of these ions over much of 
the adopted range in U, due to the large recombination rates. Nevertheless, An = models 
predict larger columns of these ions for Log U ;^ 0.7, as a result of the faster recombination 
of the L-shell ions. 

In order to illustrate the effects of uncertainties in the estimated DR rates, we 
introduced a "noise factor" by using a random number generator with a Gaussian 
distribution, with a standard deviation of 0.3 dex, suggested by the largest error in the 
polynomial fits (see Figure 1). Furthermore, the DR rates computed by Gu (2003) for 
L-shell Ar vary by ~ ± 0.3 dex over the range in electron temperatures predicted by our 
models. Each time Cloudy was run with this option, the recombination coefficient was 
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multiplied by the random number. The scale factor was kept constant for a particular 
calculation, but different random numbers were used for different calculations. We 
generated a set of 20 models at each value of U, and determined the standard deviation 
of predicted column densities. In Figure 5 we show the scatter for Fe X at three values of 
U. Clearly, uncertainties in the rate coefficients have an measurable effect on the predicted 
ionic columns. Nevertheless, the trend to peak abundance at higher ionization parameter is 
still obvious, as demonstrated by the overplot of the Fe X columns from the models without 
the estimated DR. 



3. Conclusions 

The results of this simple test indicate that the model predictions of the Fe ionization 
balance depend critically on accurate An = DR rate coefficients. Photoionization models 
of X-ray absorbers that match spectral features from abundant 2nd row elements such 
as C, N, O, and Ne will overpredict the average Fe ionization state if this process is not 
included. This can be easily demonstrated using Figures 1 - 3, by selecting an ionization 
parameter to give the proper ratio of the O VI, O VII, and O VIII column densities and 
reading off the relative Fe ionic columns for the models with and without the estimated 
DR. The apparent discrepancy in the ionic states detected in the Fe M-shell UTA in NGC 
3783 (Netzer et al. 2003), compared to those predicted by the absorber models, is likely due 
to the absence of An = rates. To illustrate this point, we regenerated the so-called low 
ionization model used by Netzer et al. to fit the soft X-ray absorption and inner shell Si and 
S lines detected in a 900 ks merged Chandra/ High Energy Transmission Grating spectrum 
of NGC 3783 (see, also, Kaspi et al. 2002). Following Netzer et al, we assumed the SED 
described in Kaspi et al. (2001) and set log U = —0.26 and N^ = 10^^-^ cm~^ (note that our 
conversion from the X-ray ionization parameter used by Netzer et al. to the form of U used 
in CLOUDY resulted in a model prediction of slightly higher relative ionization of Si and 
S than they predicted, which may be due to slight differences in the model SED; in order 
to reproduce the same columns of Si and S we adjusted the value of U downward by 10%). 
The predicted ionic columns for Fe Vlll through Fe XII, without An=0 DR, are given in 
Table 2, and the combined Fe XI and Fe XII column is ~ 3.5 greater than the combination 
of Fe VIII and Fe IX, while the peak ionization state is Fe XI. Since the strongest UTA 
features for Fe XI and Fe XII he at 16.16A and 15.98A , respectively (Behar et al. 2001), 
the UTA would peak near 16. lA , as described by Netzer et al.. Also shown are the 
predictions for the same model with An=0 DR included. In this case, the peak ionization 
state is Fe X, and the total Fe Vlll and Fe IX column are approximately equal to those 
of Fe XI and Fe XII combined. Since Fe VIII and FelX produce their strongest features 
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at 16.89A and 16.5lA, respectively (Behar et al.), and the profile would be dominated by 
Fe X (16.34 A), the predicted UTA wiU peak near 16. 4 A , as observed. It should also be 
noted that including the new L-shell rates for 3rd row elements such as S, Si, and Mg will 
alter the predicted ionization balance for those elements as well. As noted by Behar & 
Netzer (2002), inner shell absorption lines from these ions can be used as effective probes 
of high-column absorbers. However care must be taken in using photoionization models to 
determine ionization states and total column densities of the absorbers using these lines as 
constraints, at least until a full set of accurate L-shell DR rates have become available. 

Scyfcrt galaxies arc also strong soft X-ray emission-line sources (e.g., Sako ct al. 
2000; Kinkhabwala et al. 2002; Turner et al. 2003). It will be important to understand 
the connection between this component and forbidden "coronal" emission lines such as 
[Fe X] A6374, [Fe XI] A7892, [Fe XIV] A5303 and [S XII] A7611 (e.g. Kraemer & Crenshaw 
2000), especially since the optical/UV spectral images obtained with the Hubble Space 
Telescope/Space Telescope Imaging Spectrograph provide tighter constraints on the spatial 
extent and kinematic properties of the emission-line gas. Accurate DR rates are crucial for 
modeling this component and exploring the link between the X-ray and optical emission 
lines. 

While the models generated for this Letter clearly illustrate the magnitude of the 
problem, we are left with the question on how best to proceed in the absence of a full 
set of An = DR rate coefficients. Our Monte Carlo calculations (see Figure 5) provide 
an estimate of the error in the predicted column densities due to the uncertainties in the 
rate coefficients. In fact, our results suggest that comparisons between calculated and 
observed column densities be given error bars ~ 0.3 dex, assuming our rates are included 
in the calculations, due to the uncertainties in the rates and the unmodeled temperature 
dependence (as suggested by the variations in the published L-shell rates over the range in 
temperature predicted by our models, e.g., 1.8 x lO^K < T < 4.5 x 10^ K). Clearly, using 
rates estimated from other ions is far better than assuming the process docs not occur, 
which is what is done when it is neglected. Hence, we suggest that it is best to adopt these 
estimated rate coefficients with this stated error. Ultimately, new atomic data, in general, 
and DR rates, in particular, are needed to reliably compare the ionization of third and 
fourth row elements. 

S.B.K. and J.R.G. acknowledge support from NASA grant NAG5-13109. S.B.K. thanks 
Hagai Netzer for useful discussions concerning the Fe UTA lines. We thank an anonymous 
referee for useful comments. 
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Fig. 1. — The An = DR rate coefficients for L-sfiell Fe ions (crosses), L-sfiell Ar ions (stars), 
and tfie mean rates for tlie 3rd and 4tli ionization states of tfie CNO ions (diamonds) are 
sliown as a function of tlie residual cliarge of tiie ion before recombination. Tlie L-sliell 
rates are tliose from Gu (2003), for a temperature of 4.0 eV. Tlie dashed-lines are 2"°'-order 
polynomial fits for the three shells (N-, M-, and L-), which we used as the functional form 
for our estimated Fe rate coefficients. The dotted line is a linear fit to the rates, which we 
used as the functional form for the other 3rd and 4th row elements. 

Fig. 2. — Predicted ionic column densities for Fe VIII (dotted line), Fe IX (dashed line) and 
Fe X (dash-dotted line) are shown as a function of ionization parameter (U). The blue-lines 
and red-lines are the predictions with and without the estimated An = rates, respectively. 
The solid black lines are the predicted column densities for O VI, O VII, and O VIII, as 
indicated. For both sets of models, we assumed solar abundances and a hydrogen column 
density — 10^^ cm~^. 

Fig. 3. — Predicted ionic column densities for Fe XI (dotted line), Fe XII (dashed line) and 
Fe XIII (dash-dotted line) are shown as function of U. All other symbols are as defined for 
Figure 3. 

Fig. 4. — Predicted ionic column densities for Fe XIV (dotted line), Fe XV (dashed line) 
and Fe XVI (dash-dotted line) are shown as function of U. All other symbols are as defined 
for Figure 3. 

Fig. 5. — The predicted Fe X column density is shown as a function of U. The solid-hne and 
dashed-lines are the predictions with and without the estimated An = rates, respectively. 
The error bars are included to show the effect of uncertainties in the DR rates, which we 
modeled by introducing a random "noise factor" (see text). 
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Table 1. Polynomial Fit Coefficients for the Estimated Fe DR rates 



Z CI C2 C3 



17- 23 -52.5073 5.19385 -0.126099 

9 - 15 -10.9679 1.66397 -0.060596 
< 8 -3.95599 1.61884 -0.194500 
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Table 2. Predicted Fe M-Shell Column Densities (cm-^) for NGC 3783' 



Ion log(Nio„), no An = DR log(Nio„), with A N= DR 



Fe XIII 16.62 15.90 

Fc XII 16.83 16.48 

Fo XI 16.88 16.88 

Fe X 16.78 17.02 

Fe IX 16.52 16.90 

Fe VIII 15.81 16.18 



'^Predictions are for the low-ionizatiori component of the X-ray 
absorber described in Netzer et al. (2003), assuming: log(U) = -0.26, 
log(Nff) = 21.90. 



